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Abstract 


Within the factorization-assisted topological-amplitude approach, we studied the 33 charmless B(s) — 
VV decays, where V stands for a light vector meson. According to the flavor flows, the amplitude of each 
process can be decomposed into 8 different topologies. In contrast to the conventional flavor diagrammatic 
approach, we further factorize each topological amplitude into decay constant, form factors and unknown 
universal parameters. By XY? fitting 46 experimental observables, we extracted 10 theoretical parameters 
with x? per degree of freedom around 2. Using the fitted parameters, we calculated the branching fractions, 
polarization fractions, CP asymmetries and relative phases between polarization amplitudes of each decay 
mode. The decay channels dominated by tree diagram have large branching fractions and large longitudinal 
polarization fraction. The branching fractions and longitudinal polarization fractions of color-suppressed 
decays become smaller. Current experimental data of large transverse polarization fractions in the penguin 
dominant decay channels can be explained by only one transverse amplitude of penguin annihilation dia- 
gram. Our predictions of those not yet measured channels can be tested in the ongoing LHCb experiment 


and the Belle-II experiment in future. 
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1 Introduction 


Charmless hadronic B-meson decays have been of significant interest, as they can provide us an abundant source 
of information on the flavor physics, within and beyond the standard model (SM). After the very successful 
first generation of B-factory experiments, BaBar and Belle [1], the interest in this field is reinforced by the 
LHCb experiment [2] and the upcoming start of Belle-II experiment [3, 4]. In the long term run, the LHCb 
upgrade plan promise excellent future opportunity[5, 6]. The FCC-ee, as well as CEPC, the proposal for future 
electron-positron collider will give further chance for the flavor physics study [7]. In the theoretical side, 
beyond the naive factorization approach [8], three major QCD-inspired approached had been proposed to deal 
with charmless nonleptonic B decays, based on the effective theories, namely, the QCD factorization (QCDF) 
[9], perturbative QCD (PQCD) [10], and soft-collinear effective theory (SCET) [11]. The difference between 
them is only on the treatment of dynamical degrees of freedom at different mass scales, namely the power 
counting. Within these approaches, most decay processes have been studied, including the branching fractions 
and the CP asymmetries. However, the factorization for hadronic matrix elements is only proved at the leading 
order in 1/mp, with mp denoting the b quark mass. The precision of these approaches is limited to the leading 
power calculations. 

In contrast to the above approaches based on the perturbative QCD, another idea based on the topologi- 
cal diagrams and flavor SU(3) symmetry was also proposed [12], where the nonperturbative parameters are 
extracted directly from experimental data. Therefore, the extracted parameters include the effects of strong 
interactions to all orders, as well as long-distance rescattering. This idea has been used to analyze hadronic B 
meson decays extensively [13], as well as D meson decays [14]. Although no direct power expansion is needed 
in this approach, flavor SU(3) symmetry is required to reduce the number of free parameters to be fitted from ex- 
periments. As the experimental precision is better and better, the limitation of theoretical precision is retarded. 
Recently, the improved version, the so-called factorization-assisted topological-amplitude approach (FAT) [15] 
was proposed, in order to deal with the SU(3) breaking effects. By using some of the well defined factorization 
formulas to include most of the SU(3) breaking effects, the theoretical results of the two body non-leptonic D 
decays accommodate experimental data very well. Recently, the FAT approach has been utilized to study the 
two-body charmed nonleptonic B mesons decays [16]. Within 4 universal non-perturbative parameters fitted 
from 31 experimental observations, 120 charmed B decay modes were calculated. Both branching fractions 
and CP asymmetry parameters are in agreement with experimental data well. Very recently, the charmless 
Bis) — PP and B(,, — PV processes are also studied using this approach [17]. The long-standing B > nog 


and B — Kz CP puzzles can be explained simultaneously. 


In contrast to B — PP and PV decays, charmless B(;) — VV decays are much more complicated, be- 
cause more helicity amplitudes will be considered. Due to angular momentum conservation, there are three 
independent configurations of the final-state spin vectors: a longitudinal component where both resonances 
are polarized in their direction of motion, and two transverse components with perpendicular and transverse 
polarizations. For the V — A coupling of the SM, a specific pattern of the three helicity amplitudes is naively 
expected [18], such that the longitudinal polarization fraction fz should be close to unity, while the transver- 
sal contributions are suppressed by Aocp/mg. In 2004, large transverse polarization fractions (around 50%) 
of B — $K* have been measured in the experiments. Later on, some other penguin-dominated strangeness- 
changing decays, such as B — p K* and B, — $6, have also been found with large transverse polarization 
fractions. These large unexpected transverse polarization fractions have attracted much theoretical attention 
with several explanations based on the QCDF [19, 20, 21, 22, 23, 24, 25], PQCD [26, 27, 28], even on the new 
physics scenarios [29, 30, 31, 32, 33, 34, 35, 36, 37, 38]. These decays have rich observables, some of which 
are regarded as good places for testing the SM and searching for possible effects of new physics beyond the 
SM. 

To our knowledge, these decays with two vector meson final states have not been studied in the flavor 
diagram approach. In this work, we shall explore the charmless two-body non-leptonic B(s) + VV decays, in the 
newly established FAT approach. The branching fractions, CP asymmetries, as well as the angular distributions 
will be investigated. The organization of the paper is as follows: In Section 2, we give the definitions for helicity 
amplitudes, angular variables and polarization observables. The calculation of the B — VV decay amplitudes 
in FAT framework is briefly reviewed. Section 3 provides the numerical results and the phenomenological 


discussions. We will summarize work in Section 4. 


2 Framework 


We consider a B meson with four-momentum pg decaying into two vector mesons V; (m1, p1, 1]*), Vo(mo, pa, €"), 
where ņ* and &* are the polarization vectors of each final state meson. The decay amplitude can thus be de- 


composed into three parts, 
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With definite helicity, we can write the amplitudes as : 


P= af (B isti) pas) (zh (s - 3). Q) 
A~ = ef (B —> Vi(pi,ni)Vo(pa,€1)) =i (S1 FS3). (3) 


In the naive factorization, the helictiy amplitudes af -vıy are proportional to 


af = Vex (Vel Bas)v-alB) V8 (4a). (4) 
If we ignore the m? terms, the above functions can be simplified as 
A = i SE Veru fim (0), (5) 
at = iE Vem fata [o (mg +m) APY (0) + (me =m) vq). © 
A = iE Veru fam {- (mg + m)AT^ (0) — (mg — m) V" (0)) , (7) 


where An (0), api (0) and V7" (0) are transition form factors, the definitions of which are given in ref.[8]. It 


is apparent that the transverse amplitudes .c/^ are suppressed by a factor m2/mg relative to °. Due to the 


cancelation between the axial vector form factor A; and the vector form factor V, we can arrive the relations as 
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Alternatively, we can also adopt the transversity convention to get: 
AT + A AS — Gf 
b= f, ERR E E OF (9) 
v2 v2 
Then, we have 
Gr 2 ABVi 
ep = Fg Vem fami (0), (10) 
sf =  —iGrVcku fomo (mg +m) AT" (0), (11) 
sf, = iGrVcky famo (mg — mi) V^ (0). (12) 


Any of the two vector mesons in the final sates will decay via strong interaction. They form two decay planes 
with various decay angles. Thus, for any given B — VV decay, one can define five typical observables corre- 
sponding to the branching fraction, two out of the three polarization fractions fL, fj, fL, and two relative phases 


|, ox: 
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Apparently in the naive factorization, fz is expected to be close to unity from eq.(8). 


CP symmetry is violated in weak interactions, thus one expects to have different numbers for observables 
of B — VV decay with those of its CP-conjugate B decay. The CP averaged decay rate and the CP asymmetry 


are then defined as 
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where p is the 3-momentum of the final state vector meson in the rest frame of the B meson. Observables id 


and oF are also defined as in (13), and CP averages and asymmetries are calculated by 


f= 0M. as fe — fr (16) 
i ES 1. 
(h = L, ||, L) for the polarization fractions and 
CEPMCET MEE =5(0F 08), (m 


for the phase observables $; and A@,. Unless otherwise indicated, for each observable quoted we imply the 
average of a process and its CP-conjugate one. 
For the neutral B meson decays, if f = f, the time-dependent of CP violation is defined through 


ree) > f) -r*(8*(r) 2 f) 
TB’ (t) ^ f) - T^(B9(r) ^ f) 


A^ (t) — 


cp 


E Ssin(Ampt) — C/cos(Ampt), (18) 


where Amg > 0 is the mass difference of the two neutral B meson mass eigenstates. The definitions of two 


parameters C f and si are given as: 


hy2 l 
= Tem ^ = Im (19) 
f qM 1-- AP? 
where the parameters standing for the mixing of neural mesons read 
Aj VV 
ARBí I X Bpod. (20) 


OPAR p VeV 
Al is the decay amplitude of B^ — f" and A^ is the the CP-conjugate one. Obviously, the CP violations in 
B — VV are much complicated than those of B — PP and B — PV modes. 

Now, we shall introduce the FAT approach. All of the hadronic B decays are weak decays, which are 
perturbatively calculable. However, all the initial and final states are hadrons, which involve non-perturbative 
QCD effects. The factorization of perturbative and non-perturbative QCD is limited to certain power expansion 
of 1/mp, which give limited accuracy of theoretical precision. The topological diagram approach does not rely 
on the QCD factorization, but group all the decay amplitudes by different Feynman diagrams according weak 
interaction. That means only factorization of weak interaction from strong interaction is required. The QCD 
corrections to each weak decay diagram are extracted from experimental data, including all perturbative or 
non-perturbative ones. In this case, this approach is a kind of model independent method to deal with hadronic 
B decays. Among these weak Feynman diagrams, we have four kinds belong to the process induced by tree 


diagram, which should be the leading contribution shown in Fig.1. They are denoted as 


d,s 


gi 


(o) E (a) A 


Figure 1: Topological tree diagrams contributing to charmless B — VV decays: (a) the color-favored tree 
emission diagram, T; (b) the color-suppressed tree emission diagram, C; (c) the W-exchange diagram, E and 


jm (d) the W-annihilation diagram, A. 


T , denoting the color-favored tree diagram with external W emission; 


C, denoting the color-suppressed tree diagram with internal W emission; 


E, denoting the W-exchange diagram; 


A, denoting the annihilation diagram. 


In the conventional diagrammatic approach, the numerical results of each topology diagram can be fitted 
directly from the experimental data, by using the SU(3) symmetry. It is well known that the breaking of 
SU(3) symmetry can reach 20 — 30%, which indicates that the prediction power of conventional diagrammatic 
approach is limited. For the decay processes dominated by the T-type diagram, such as B + D7, the amplitudes 
can be expressed by the products of transition form factor, decay constant of the emitted meson and the short- 
distance dynamics Wilson coefficients [16], which has been proved in QCDF, PQCD, and SCET. Results of all 


the three approaches agree with each other and with the experimental data well. Similarly, in our B(j + VV 


decays, the T amplitudes of all three polarizations are expressed by 


T? = i Meran (uU) fompAe™! (0), (21) 
T! = —iGr Verma: (H) fomo (mg +m) AP" (0), Q2) 
T+ = iGgVckyai (LU) fzm (mg — mi) V^ (0). Q3) 


In contrast to the conventional diagrammatic approach, where decay amplitude needs to be fitted from experi- 
ments, it is obvious that no free parameter is introduced in the T amplitude. a (u) = Co(u) + Ci(u)/3 is the 
combination of the effective Wilson coefficients of four-quark operators [39], where u = mp/2 = 2.1GeV is 
the factorization scale. At this scale, a1 (2.1GeV) = 1.05. The SU(3) breaking effects are automatically kept, 
by different decay constants and form factors for different processes. 

For the color-suppressed tree diagram shown in Fig.1(b), the nonfactorizable contribution is dominant, thus 
no factorization formula is given. The decay amplitude and strong phase need to be fitted from experimental 
data. However, we can still factorize out the corresponding decay constant and form factor to keep the SU(3) 


breaking effects as an approximation: 


Gp " 
C? =i V, 0 itc fym%A™ (0 24 

va VexmXc fompAg '(0) (24) 
cl = —iGrVcgm%le e'* fymy (mg +m Ae" (0) (25) 
C. =iGrVcxuxee!% fom (mg — mi) V9" (0). (26) 


Similarly, for the W-exchange diagram shown in Fig.1(c), we factorize out the meson decay constants to 


characterize the SU(3) breaking effects. The decay amplitudes for the three polarizations are then written as: 


_Gr hO 1 

E? — i— Vcky fa fi fom xpe E — (27) 
v2 fo 
I1 

El = —iGrVckw fa fi fomaxne » f Q8) 
i 

L = iGFV, BHE eE 7 2 
= iGrVcxu fa fi fom xg e A (29) 


where the amplitude yg is dimensionless with normalization to fo: As for the W-annihilation diagram A shown 
in Fig.1(d), its contribution is too small to be considered in the present work. In practice, even if we keep this 
contribution in our fitting program, we cannot get a stable solution for them, since the experimental precision 
is not good enough. 

There are also QCD-penguin diagrams shown in Fig.2. Although they are loop diagrams, which should be 
suppressed comparing with the tree level diagrams, they are enhanced by the large CKM matrix elements and 


large top quark mass. There are also four types of them: 


(c) Pr (d) Pa 


Figure 2: Topological penguin diagrams contributing to charmless B — VV decays: (a) the color-favored QCD- 
penguin diagram, P; (b) the flavor-singlet QCD-penguin diagram, S; (c) the time-like penguin annihilation 


diagram, Pg and (d) the space-like penguin annihilation diagram, P4. 


P, denoting the QCD penguin diagram; 


S, denoting the flavor-singlet QCD-penguin diagram; 


e Pr, denoting the time-like penguin annihilation diagram; 


P4, denoting the space-like penguin annihilation diagram. 


Similar to the color favored tree diagram T, the QCD penguin diagram P has the same V — A structure, 
which can be proved factorization for all orders of Œs expansion. Therefore, we do not need to introduce free 
parameter for them. The decay amplitudes are again expressed by the products of transition form factor, decay 


constant of the emitted meson and the short-distance Wilson coefficients: 


p — MA fom A8“ (0), (30) 
Pl = —iGrVcgma4( U) fom (mg 4- m1) AP^ (0, ) (31) 
PŁ =iGrVcxmaa(LL) fomo (mp — m1) V9" (0). (32) 


Similar to the color suppressed tree diagram C, the flavor-singlet QCD penguin diagram S are non-factorizable 


that expressed by 
Gr “40 
S9 = i — Vcka xL e*s f,m2AB" (0), 33 
aoe fomgAg | (0) (33) 
sl = —iGrVcru%le e fam; (mg +m) A^ (0), (34) 
SŁ = iGrVexm xi e fom (mg = mı) yen (0). (35) 


Similar to the W-annihilation type diagram, the contribution of time-like penguin diagram shown in Fig.2(c) 
is also negligible, which can not be fitted from the experimental data at the current precision. Lastly, the penguin 


annihilation diagram (space-like) contributions P4 are expressed as: 


G ip | 
P} = i Vex fol fanis e a SR 
p! = —iGrVcru fsfi fama X), e e 3 fg Gn 
p 
1 
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Now, we have 12 magnitudes and 12 phases to be fitted simultaneously from the the experimental data. 
With the fitted parameters, we can predict all other B; > VV (q = u,d,s) decays with branching fractions, 


polarization fractions, and CP asymmetry paramters. 


3 j?Fitting and Numerical Results 


To characterize the flavor SU(3) breaking effects in our calculation, we need input parameters of various meson 
masses [41] and decay constants that are summarized in Table 1. Unlike the meson masses, the value of decay 
constants is not known in experiment, but can be given from theoretical calculation, such as QCD sum rules 
[46], Bethe-Salpeter equation [47], and lattice QCD [48], and we taken the value from [17] with 5% uncertainty. 

For the calculation of color favored tree diagram T and QCD penguin digram P, we also need the input of 
form factors for B — V transitions. There are many calculations for form factors, such as light-cone sum rules 
[49], perturbative QCD approach [50], lattice QCD [51] etc. The central values we used in this work of the 
transition form factors at q? = 0 are shown in Table 2. To estimate the theoretical uncertainty of the numerical 
results in our calculation, we include the uncertainties of all form factors as large as 10%. In fact, they are one 
of the major sources of theoretical uncertainty in our numerical results, especially for processes dominated by 
the color favored tree (T) diagram. Since the final state meson mass is small comparing with the large B meson 
mass, the q? dependence of form factors will be neglected. In fact, the effects of g? dependence to numerical 


results are negligible, which has been indicated in B — PP, PV decays [17]. 


Table 1: The mass and decay constant of meson (in units of GeV) 


Meson Mass Decay Constant 
pe 5.28 0.190 
B? 5.36 0.225 
p 0.77 0.213 
0.78 0.192 
Q 1.01 0.225 
K* 0.89 0.220 


Table 2: The transition form factor of B — V at q? = 0 


Bop B—oK'* B—5o B,—cK' B;>ọ 


V(0 033 041 029 0.31 0.42 
(0) 032 038 0.28 0.36 0.44 
(0) 025 0.29 0.22 0.23 0.31 


A1 


For the CKM matrix elements, we adopt the Wolfenstein parametrization, which are given as [41] 


A — 0.22537 +0.00061, A = 0.814* 0024, p — 0.117 2:0.021, 7 = 0.353 0.013. (39) 
The life time of B mesons are taken from the particle data group [41], given as 


Tp- = 1.641 ps, tgo = 1.519ps, Tgo = 1.497ps. (40) 


3.1 The Y? fit for theoretical parameters 


Unlike B > D“)M and B > PP(V) decay modes, the amplitudes of B — VV modes are much complicated, 
because each decay process has three polarization contributions, which means that the number of parameters 
will increase threefold. From previous section, we do not introduce any new parameter in color favored tree 
diagram T and QCD penguin diagram P. For the color-suppressed tree diagram C, the flavor-singlet QCD- 
penguin diagram S, the W-exchange diagram and QCD-penguin annihilation diagram P4, we have all together 


OIL OG), 0L olL OIL LOL 0,1) 0, ||, 
24 parameters xe. ; el , xel ; el E ; orl , xe and op! . So many free parameters are 


difficult to be determined from the limited number of experimental measurements. It will also decrease the 
prediction power of this FAT approach. As indicated from the QCD factorization approach [24, 25] and the 


perturbative QCD approach [28] calculations, the color suppressed tree diagram C, W exchange diagram E 
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and the flavor singlet QCD penguin diagram S are dominated by longitudinal polarization contributions. For 


= xt = 0. Therefore, the only transverse polarization amplitudes to be fitted are from 


simplicity, we here drop the negligible transverse contributions of these topological diagrams to set y, 
0, x! ess and x 
the penguin annihilation diagram. According to the power counting of this diagram [20], the negative helicity 
amplitude is Chirally enhanced while the positive helicity amplitude is still suppressed (to be neglected). This 
results in a relation of XP, RI Xe, . Thus, there are only 10 universal parameters left, which will be fitted by 
experimental data. 

In the experimental sides, after the first decay modes B — @K* were measured, more and more observables 
have been measured, involving the branching fractions, CP asymmetries, polarization fractions, and relative 
phases of helicity amplitudes. Since some observables are measured with very poor precision, the data with 
less than 3o significance will not be used in our fitting program. Then, we have 46 experimental data, involving 


18 branching fractions, 20 polarization fractions, 6 relative phases, and 2 direct CP asymmetries. 


In the Y? fitting, in term of 46 experimental observables x; + Ax; and the corresponding theoretical predic- 


tions x^, the y? function can be defined as 


46 / th 2 
2. Xj —i 
e=¥ (HS). an 


With the amplitudes and data, we can extract the 10 parameters by minimizing the x?. The best-fitted values of 


the parameters are given as: 


Xe = 0.23 + 0.05, oe = 0.48 +0.29; 72 = 0.082 +0.026, ó = 1.69 +0.16; 


Xs = 0.018 +0.003, 42 = 1.29+0.22; xP =0.012+0.002, op, = —0.07 +0.18; 


xi = 0.0098 + 0.0003, ¢}* = —0.21 +0.09; (42) 


The %?/d.o.f = 82.0/(46 — 10) is 2.28. Compared with the corresponding fitted values of B — PP and PV 


decays [17], the magnitudes of the longitudinal polarization X? (i = C, E,S, P4) are at the similar size. 


3.0 Numerical Results and Discussions 


For the convenience of later discussion, we shall class 33 decay modes into 5 categories. The first category is 


T -dominated decays, involving four decay channels, B 4 p'p .B —p*p9,B- — pto and B 2p K*. 
Five decay process B —> p?p’, B —> p?o, B — p°K* and B — wK*° dominated by C diagram, fall into the 
second class. The third class is dominated by QCD penguin diagram P, with eleven decays B —> pK*, @K*, 
oK*, K*K* and B, => K*K*, @K* and ġġ. There are three more decay channels B^ — K*°K*~, B > KOK” 


and B — @K*°, which are also P dominated. But their branching ratios are highly suppressed comparing with 
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the previous 11 channels, since they are mediated by b — d transition instead of b — s transition, which are 
CKM suppressed. The fourth class includes B- > p- $, B? > p°¢, B? > c, B, > pọ and B, > «6, to 
which the flavor-singlet QCD-penguin mainly contribute. The last group is five decay channels dominated by 
annihilation type diagrams: B? + K*^K*^, B, > p* p^, B, > p?p’, B, — p°@ and B, > c. 

In Table 3, we list the branching fractions, and the theoretical errors correspond to the uncertainties due 
to variation of (i) the fitted universal y values, (ii) the heavy-to-light form factors and (iii) the uncertainty 
of decay constants. The error of the variation of the CKM matrix elements is negligible. We note that for 
other observable, we have combined these uncertainties by adding them in quadrature and show the resulting 
uncertainty, due to to the space limitations in the tables. The decay B — $$ is absent in our tables, because it is 
a pure annihilation type process induced only by time-like penguin diagram Pg, which contribution is neglected 
in this work, though its branching fraction is estimated to be at the order of 1078 based on PQCD approach 
[42]. Comparing our predictions with experimental data, one can find that our results can accommodate the 
data well, within uncertainties from both theoretical and experimental sides. For those decays that have not 
been measured, our prediction can be tested in the ongoing LHCb experiment and the forthcoming Belle-II 
experiment. 

For the color-suppressed decay mode B 5 p9 o, its branching fraction is predicted to be 0.08 x 1076 and 
0.4 x 1076 in QCDF [23] and in PQCD [28], respectively. However, in this work, we predicted its branching 
ratio to be 1.48 x 1076, which is about 18 times larger than the result of QCDF. Moreover, in both QCDF and 
PQCD, the longitudinal fraction is calculated to be about 67%, which is smaller than our prediction (87%). 
The measurement of this mode in future will help us to understand the dynamics of this decay mode. For 
B — œo, although the experimental data agree quite well with our value, but there exist larger uncertainties in 
both experimental and theoretical results. So, the precise measurement of B — 00 is necessary, too. 

The polarization fractions, as well as relative phases, are shown in Table 4. One can find that for the first 
four tree-dominated processes, they fully respect the helicity hierarchy that are dominated by longitudinal po- 
larization. The major theoretical uncertainties here are from the heavy-light form factors and the CKM matrix 
element |V,,|. On the contrary, for the color suppressed tree diagram (C) dominant decays, their branching 
fractions become much smaller due to the absence of T-type diagrams. Correspondingly, the largest uncertain- 
ties are from the xc, in particular from the transverse polarizations, which have also been confirmed in QCD 


factorization [21, 23]. For the three B — pp decay modes, the decay amplitudes can be read: 


V2A(B- > pp?) 2 T +C; (43) 
A(B’ 2 p*p~) -T -E 4 P Pi; (44) 
A(B° > p?p*) = -C+E +P +P. (45) 
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Table 3: Branching fractions and the direct CP asymmetries of B + VV decay modes. The experimental data 


[41] are also given for comparison. 


Branching Fraction / 1076 


Acp / percent 


Cras eas Mad Theory Exp Theory Exp 
B pp | 21.741.8+4.24+2.2 24.0419 |0 —5+5 
Bp pt | 29.541.945.443.0 28.3421 | —8.10 +2.94 
B-—5p-o | 18.241.54+2.8+1.6 15.942.1 | —3.4545.38 | 2049 
B? = K*+p- | 38.640.1+7.343.9 —10.9+3.0 
B > p?p? 0.94+0.46+0.11+0.14 | 0.9740.24 | 49.7 +13.4 
B" > p?o 1.48 +0.71 +0.06+0.20 | «1.6 —38.5 +13.6 
C |B > oo 1.20+0.49+0.12+0.18 | 1.2404 | 28.1413.8 
B? — K'9p? | 1.1840.39+0.2140.12 | <767 4.94183 
B? Ko | 0.9740.33+0.16+0.10 32.2+16.0 
B- >p- K” | 104+1.6+1.7+1.1 924-15 |1.004017 | -14£16 
B- —p?K*- | 5.83+40.66+0.7640.65 | 4.6+1.1 34.6 8.3 3113 
B — pK | 5.0940.7540.82+0.53 |3.9+1.3 |-—0.6+4.0 | -6+ 
B? — ptK* | 10.541.341.4412 10.342.6 | 343463 | 21415 
B- + @K*- | 4.2440.7040.32+40.51 | «74 —30.1 +13.8 | 29435 
P |B oK | 3.1040.7540.2740.38 |2.0+0.5 | -11.744.0 | 45425 
B- + @K*- | 9.3141.9041.8340.97 | 1042 1.0040.27 | -1+48 
BOK”  |864-1.76:1.70-0.90 | 1040.5 | 1.0040.27 +4 
B SKK"? | 14.942.041.9423 2847 0.78 +0.19 
B? > K*-K*+ | 15.941.74£1.742.6 21.147.1 
B. — 09 26.444.844.543.8 19.343.1 | 0.83 +0.28 
B-  K'K*- | 0.66£0.1040.1340.08 |1.2+05 | —24.8+2.6 
P |B? KVK” | 0.6140.09+0.12+0.07 —24.84-2.6 
B? — @K* | 0.70+0.11+0.13+0.08 | 1.1340.30 | -17.345.6 
B-—óp- | 0.0640.0240.0140.01 | <3.0 0 
B^ > op? 0.0340.0140.0140.00 | <0.33 0 
S |B 560 0.02 +0.01 +0.00 +0.002 | < 0.7 0 
B? > op° 0.07+0.03+0.0140.01 | <617 0 
Bo > ġo 3.69+1.19+0.7440.37 —15.0+7.0 
B" — K*+K*- | 1.4340.9140+0.29 «90 0 
B'—p-p* | 0.1040.0640+0.02 0 
E(Pg) | Bo p?p?  |0.05::0.08::040.01 < 320 0 
B^ 5 pw 0.08 £0.05 +0 +0.01 0 
B > oo 0.03 2-0.02 3-0 2- 0.01 0 
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Table 4: The polarization fractions and relative phases of B —> VV decay modes. The experimental data are 


taken from [41]. 


ÍL/ percent f1/ percent $j/ rad $,/rad 

Decay Mode 

Theory Exp Theory Exp Theory Exp Theory Exp 
B- + p7p’ 95.5+1.1 | 9541.6 | 2.22+0.64 —0.09 + 0.05 —0.09 + 0.05 
B > p-pt 92.6--1.6 | 98.8+2.6 | 3.65+0.91 —0.27 + 0.08 —0.27 4: 0.08 
B >p o 92.741.4 | 906 3.60+0.76 —0.23 + 0.07 —0.23 +0.07 
B= K*tp- | 944412 2.74 3: 0.64 —0.08 + 0.03 —0.08 + 0.03 
B = p?p? 81.7410.8 | 60423 9.21 3: 5.50 —0.04+ 0.44 —0.03 + 0.44 
B^ = plo 82.7 9.5 8.68 +4.82 0.98 +0.22 0.98 +0.22 
B > oo 92.2+3.6 3.94 + 1.85 —1.46+0.26 —1.45+0.26 
Be — Kp? 79.8 +8.0 10.2+4.1 —0.94 + 0.28 —0.94 + 0.28 
B? 5 Kg 71.9 9.2 11.22: 4.7 —0.73 +0.31 —0.73 +0.31 
B-—p-K | 46042129 | 48+ 272 37.0 2.07 3: 0.22 2.08 + 0.22 
B-—p?K*- | 40.7410.6 | 78x12 29.843: 5.9 2.24 3: 0.20 2.24 3: 0.20 
B > ook? 48.7412.3 | 40414 25.8 +6.7 2.08 +0.21 2.09 +0.21 
B? => ptk*- | 38.9+11.3 | 38413 30.8 46.3 2.18 +0.22 2.18 +0.22 
B-— @K*— | 29.9+6.8 | 41+19 35.3445 0.02 + 0.85 0.03 + 0.85 
B 5 ek"? 29.4417.5 | 69413 | 35.694 ~2.62+£0.53 —2.61 +0.53 
B- > oK*- 48.0+16.0 | 505 25.94+8.6 | 20+5 2.4740.27 | 2.34£0.18 | 2.47+0.27 | 2.58+0.17 
B 5k? 48.0+16.0 | 49.7+1.7 | 26.048.6 | 22.441.5 | 2470.27. | 2.434+0.11 | 27:027 | 2.5340.09 
B? => KOKO | 343412.6 33.2+6.9 | 3811 2.10 +0.23 2.10 +0.23 
B? + K*-K*+ | 30.9+10.4 34.945.8 2.19+0.22 2.20 +0.22 
B 566 39.7416.0 | 3621.4 | 31.248.9 | 30.941.5 | 2.5340.28 | 2.55+0.11 | 2.5640.27 | 2.67::0.23 
B- —K"K*- | 58.3411.1 | 75421 20.8 4: 6.0 2.1023- 0.20 2.09 + 0.20 
B — KR | 58.3411.1 | 80413 20.8 + 6.0 2.10+ 0.20 2.09 + 0.20 
B? = oK” 38.9414.7 | 51417 3148.1 2.52+0.27 2.55+0.27 


The three decay amplitudes make an isospin triangle. For illustration, we list the numerical results of longitu- 


dinal polarization for each topological diagram of these decays, 


|r8?P:9:| ; |CP>PLPL] ; |EP>PLPL] ; pB>PLPr] ; |pP>PHPL| — 1 : 0.22: 0.21: 0.14 : 0.08. 


(46) 


For the decay B7 — p p, although the absolute value of the C diagram is suppressed, it can enhance the 


magnitude of the decay amplitude by 2096. With the larger contribution from C diagram, the large branching 


fraction of B^ > p?p’ is also explained. By including the transverse momenta of inner quarks, the n 


p?p’ has been calculated in perturbative QCD approach in ref.[28], where they got a rather smaller branching 


fractions (0.27 x 1076) and a smaller longitudinal fraction (12%). In fact, there is a large discrepacy between 


experimental measurements from BaBar and Belle, and the number we quoted in this work is the naive averaged 


value. So, it is very important to have a refined measurement of the branching fractions and longitudinal 
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Table 5: Amplitudes of each diagram of B — @K* (x 10~8GeV) 


P S Py 
A9 —1.36--501i 1.39—043i  —0.21— 2.29i 
Al  —0.34-122i 0 4- 0i —0.60 — 2.60i 
At —0.33+1.23i 0+0i —0.60 — 2.60: 


fractions of B — pp decays to draw the final conclusion. 

For B, decays dominated by T and C diagrams, they have the same manner as B decays, for example, 
the color-allowed decay B — p K** has large branching fraction and large longitudinal fraction, while the 
branching fraction of color suppressed B — K*9p? (m) decays is a bit smaller and the transverse polarizations 
are about 20%. Comparing the branching fractions of B — p K** with B > p p^, we find that the former 
is larger than the latter one. The reason is that the form factor Ae is larger than A by 13%. Considering 
the life time difference, the large gap between these two branching fractions can be well understood. Since 
the order of these three branching fractions is about 1076, they should be measurable in the running LHCb 
experiment. 

We now discuss the 14 decay modes dominated by QCD penguin diagrams P. Among these decays, 11 
decays induced by b — s transition have branching fractions up to 10^? due to large CKM matrix elements 
\V:pV,<|. Some of them have been measured precisely in the experiments, including branching fractions, polar- 
ization fractions, and even CP asymmetries. The first measured one, and also the most well measured channel 
is B — @K* modes that are induced by b — sss transition. In this decay, the magnitudes of QCD penguin 
diagram P, the flavor-singlet penguin diagram S and the penguin annihilation diagram Py, are at the same order 
magnitude. For illustration, numerical results of each diagram are provided in Table 5. It is easy to see that 
the penguin annihilation diagram has a very large transverse polarization contribution. In QCDF [21, 23], in 
order to increase the effects of annihilation diagrams, the free parameters p4 and $4 introduced for the power 
suppressed penguin annihilation diagram are required to be very large. In the so-called PQCD approach, the 
large effects of annihilation are arrived, by including the transverse momenta of inner quarks. So, we then 
conclude that the larger transverse polarizations in B — $K* arise from the annihilation diagrams. 

Another decay mode induced by b — sss is Bs —^ $9. From Tables. 3 and 4, one can see that although 
our estimation of branching fractions agree with data within uncertainties, our center value is a bit larger than 
the experimental data, and the predicted polarization fractions are in agreement with the data. The acceptable 
divergency in branching fraction is due to the larger form factor AU we adopted, which in fact is also related 


to the branching ratio of B, + $K*. If we consider the B, — $6 alone, the present data favor a smaller AD. 
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In this point, the precisely calculations of form factor in lattice QCD and other effective approach are needed. 

ForB —p 6, B > p°@ and B. Q9? $, they are dominated by the flavor-singlet QCD penguin diagram. 
These tree decay modes have an identical amplitude S. Compared with the longitudinal amplitudes, both 
transverse amplitudes are so small that can be neglected safely, so that the longitudinal polarizations are about 
10096. It should be emphasized that the neglected electroweak penguin will enhance the negative helicity 
amplitudes, and the longitudinal polarization fraction will decrease to 70% [44]. For decays B — p° and 
B — w°, they have large branching fractions as they are governed by the larger CKM elements. In fact, 
because we drop the electroweak contribution away, and the flavor singlet penguins is cancelled, the B — po 
decay is a mode with only color-suppressed tree diagram (C) contribution. While for the B — wo, both C 
and S contribute, so its branching fraction is larger. Honestly, for this decays, there is another larger uncertainty 
we have not included here. In this work, we have assumed the ideal mixing, which means that there is no s5 
component in the @ component. As we mentioned above, B, — $$ has a large branching fraction, even a small 
mixing angle will affect the predictions remarkably [45]. 

In Table 4, as we expected, the longitudinal polarization fractions of the tree diagram dominant decays 
are predicted to be near unity with errors in the (5 — 10)96 range. The CP asymmetries in the longitudinal 
polarizations of these decays are less than 5%, as shown in Table 6. Although the CP asymmetries of the 
perpendicular polarizations of these decays shown in Table 6 are large, they are difficult to measure, since 
their fractions are too small. For the decays controlled by the C diagram, although the longitudinal polarization 
factions become smaller, they still play the primary roles with large uncertainties. Furthermore, the uncertainties 
of Ap is large, though some of them can reach 10%. Our theoretical result of fr, for the B^ — p p? is a bit 
larger than the data, but the situation of Bs p*p is in the opposite direction. 

Because each B — VV decay has three polarizations, the possible time-dependent CP violation is compli- 
cated, itis very hard to measure them precisely in the experiments. For the Tree-dominant decays, the transverse 
parts can be neglected, and the measurement of time-dependence CP violation of these decays becomes plausi- 
ble. In this work, we calculated the C; and S; of the longitudinal parts of 8 —> p'p- and BS p9p?, as shown 
in Table 7. Obviously, our results agree with experiment, though there are large uncertainties in both theoretical 
and experimental sides. Also, we present the numerical results of B? + (K*- K**); and B? — (@@)z, which 
may be measured in future, as both modes have large branching fractions and the final states are easy to be 
identified. Noted that the precise measurement of Crs and sb p Will help us to determine the CKM angles @ and 
y [43]. 

We then come to decay modes B — p K* and B, — K*K*. For the pure-penguin process B^ — px". 


which is similar to the decays B — $K"*, the penguin annihilation diagram will give large transverse polarization 
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Table 6: The predictions of the CP asymmetries of longitudinal and perpendicular polarizations, as well as 


relative phase differences of B — VV decays. 


Decay Mode AQ pl percent Acpl percent | Aój/ rad A0 / rad 
BS pp’ [0 0 0 0 
B'—p-p* | 130+0.54 | -16.348.2 | -0.4140.05 | —0.41 +0.05 
B-+p-o | 2.384086 | —30.2411.6 | -0.7040.07 | —0.72 +0.07 
Bo K*+p- | 0.9140.45 | -15.449.5 | —0.5240.06 | —0.54+40.06 
B" > pop? 105496 | —46.9+ 13.9 | 1.89+40.19 1.89 2- 0.19 
B? > plo —8.68+7.72 | 41.6413.3 | —1.4740.09 | —1.47-£0.09 
B" > wo 2.1041.87 | —24.7414.1 | -1.4340.07 | —1.43 +0.07 
B? — Ko? | 0.4744.69 | —1.89+18.3 | 2.0340.10 2.02 +0.10 
B? >Ko =| 8.374£7.28 | -29.5416.3 | —1.58 +0.13 | —1.58 +0.13 
B-p- K” | 1.404056 | —1.19:0.19 | —0.0140.00 | —0.01 +0.00 
B- + p°K*- | 35.0419.8 | —42-:9.0 | -0.8240.15 | —0.82+0.15 
B > pK" | -0.4144.3 | 0.3944.06 | 0.13+0.04 0.13 +0.04 
B > ptK* | 37.2418.9 | -23.846.9 | —0.66:0.12 | —0.65 +0.12 
B- > oK'- | -93.4425.9 | 39.6413.0 | 2.09+0.90 2.09 +0.90 
B? > oK” | -27.7419.4]| 11.544.0 | -0.0440.12 | —0.04+0.12 
B- > 6K*- |126-0.71 | —1.16+0.30 | —0.0240.00 | —0.02 +0.00 
B+ @K  |1265071 |-1.16-0.30 | —0.0220.00 | —0.02 +0.00 
B? — KOK | 1.8140.69 | —0.94--0.20 | —0.005--0.003 | —0.004 + 0.003 
B? > K*-K*+ | 32.6420.2 | 148-73 | -0.7040.14 | —0.70+0.14 
B > ọọ 1.55+0.85 | —1.02+0.29 | —0.0140.00 | —0.01 +0.00 
B —K"K*- | -20.2+8.0 283-33 0.16 +0.04 0.16 +0.04 
B? > KK | -20.248.0 | 28.3+43.3 | 0.16+0.04 0.16 +0.04 
Bl oK*® | -32.9415.0| 21.0+5.7 | 0.27+0.08 0.27 +0.08 


Table 7: Prediction of the time-dependent CP Violation(%). 


Sp Cr 

Decay Mode 
Theory Exp Theory Exp 

B? — (p- p^) —3.67 +3.02 -6+17 6.80 +3.12 —5+13 
p? => (p9p9); 41.7 X 27.8 30+70 —57.2 +17.4 2090 
B? > (m0); 15.7 +13.3 —30.0+15.1 
B)—(K*-K**), 809-144 —50.2+ 20.7 
B9 > (00): 2.16+0.76 —2.39+0.82 
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fraction. For the B- — p°K*~, to which the tree operators also contribute, the destructive interference between 
tree and penguin operators reduce the longitudinal amplitude. So, the smaller longitudinal polarization fraction 
of B- — p?K*- is obtained. The large longitudinal polarization fraction fy, of this decay is only measured by 
BABAR experiment. We hope Belle or Belle II experiment can help to resolve this puzzle. Due to the factor 
1//2, the branching fraction of B7 — p°K*~ is about half of that of B^ — p^ K “0 The analysis and the result 
of the modes B + co K* and B, — K*K* should be similar to those of B — p K*. From Tables. 4 and 6, we find 
that for all penguin dominant decays f, ~ f, 6j + 6) and Aój ~ Ag, , which indicates that the positive-helicity 
amplitudes are about zero. In fact, due to the suppression of leading QCD penguins, pK* final states have also 
been used to prob electroweak penguin effect [21], however, this kind of contributions have been neglected in 
the present work due to not enough experimental data. 

As for the last five pure annihilation type decays, they all have two kinds of contributions: the W exchange 
diagram (E) and the time-like penguin annihilation diagram (Pr). As discussed in previous section, there are 
not enough experimental data to determine the amplitude of time-like penguin annihilation diagram (Pz). In our 
fitting, we have to set it to 0. Since all the B, decays in this category are dominated by this Pg contribution except 
B — p? o, due to the small CKM matrix elements in W exchange diagram (E), the branching fractions of these 
decay modes are not stable in Table 3 with large uncertainties. Only B 5 KK- decay has a relatively larger 
CKM matrix elements in the W exchange diagram (E), which makes a relatively larger branching fractions, 
but still with large uncertainty. And the pure W exchange diagram (E) channel B — p?o also have large 
uncertainty because of the large error of %2. Therefore, we conclude that the current experimental data can not 
help us to make predictions on this kind of decays, but waiting for the running LHCb experiment, Belle-II or 


other future colliders. 


4 Summary 


In the work, we preformed analysis of 33 charmless two-body B(s) + VV decays within the factorization- 
assisted topological-amplitude approach. In contrast to the charmless B — PP and B — PV decays, more 
parameters (triple number in principle) are needed to describe the three polarization amplitudes of B(,, > VV 
decays. However, with the current 46 experimental data, we can only fit 10 universal parameters of them. For 
the decays with large transverse polarization fractions, such as the penguin diagram contribution dominated 
decays, we need only one transverse polarization amplitude in penguin annihilation diagram to explain all the 
polarization data. We calculated many decay modes not yet measured, involving the branching fractions, the 


polarization fractions, CP violation parameters, as well as relative strong phases. These results will be tested in 


18 


the LHCb experiment and future Belle-II experiment. 
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